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Anomalous metal- insulator transition observed in filled skutterudite CeOs4Sbi2 is investi- 
gated by constructing the effective tight-binding model with the Coulomb repulsion between f 
electrons. By using the mean field approximation, magnetic susceptibilities are calculated and 
the phase diagram is obtained. When the band structure has a semimetallic character with 
small electron and hole pockets at F and H points, a spin density wave transition with the or- 
dering vector Q = (1, 0, 0) occurs due to the nesting property of the Fermi surfaces. Magnetic 
field enhances this phase in accord with the experiments. 
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Filled-skutterudite compounds RM4X12 (R: rare 
earth, M: Fe, Ru, and Os, and X: P, As, and Sb) with 
the unique body-centered cubic structures have attracted 
much interest due to wide variety of physical phenomena, 
e.g. magnetic and quadrupole ordering, Kondo insulating 
behavior, unconventional superconductivity, etc. 1,2 

Recently, anomalous phase transition from paramag- 
netic metal (PM) to spin density wave (SDW) state was 
found at T =0.8K in CeOs4Sbi2. Surprisingly, the phase 
boundary is shifted to higher temperatures by applying 
the external magnetic field. 2,3 Such increase of the tran- 
sition temperature have been often observed in the sys- 
tems with quadrupole-ordering, 4 but the magnetic mea- 
surement 5 suggests that the crystal-field (CF) ground 
state of f electron is T7 which has no quadrupole mo- 
ment, and the specific heat measurement 2 shows only 
a small entropy change at the transition, excluding lo- 
cal moment ordering but suggesting an SDW or CDW 
transitions with an opening of a gap over a small por- 
tion of the Fermi surfaces of itinerant electronic states. 
The neutron scattering experiment 6 was done but the 
CF was not clearly observed. The NQR measurements 7 
claims the first-order transition since the hysterisis was 
observed, but the specific heat data 2 clearly show peaks 
at the transition indicating the second order transition 
with some short range order. 

CeOs4Sbi2 shows semiconducting behavior in the re- 
sistivity experiment with a small charge gap (E g ~ 10K) 
at low temperatures. 5 The overall structure of the op- 
tical conductivity and its temperature dependence show 
similar behaviors to those of CeRu4Sbi2 which exhibits 
a pseudo-gap and the mid-infrared peak with strong 
temperature-dependence. 8 ' 9 Although the band struc- 
ture of CeOs4Sbi2 obtained by the band calculation is 
quite similar to that of CeRu 4 Sbi2, 10 CeOs4Sbi2 shows 
semiconducting behavior, whereas CeRu 4 Sbi2 is metallic 
in experiment. 5 

In this letter, we construct the effective model with 
the characteristic feature of the filled skutterudite com- 
pound CeOs4Sbi2 and investigate the anomalous metal- 
insulator transition observed in this material. 

The basic crystal structure of filled skutterudites is 
that the X12 clusters, filled by R in the center, form the 



body-centered cubic (BCC) lattice. The f orbitals of rare- 
earth atoms strongly hybridize with the p orbitals in X12 
clusters close the Fermi level. The X12 cluster has the Th 
symmetry, but its effect on J— 5/2 f 1 state in Ce is the 
same as Oh- 

In order to investigate the properties of Ce- 
skutterudites, we previously proposed a simplified tight- 
binding band model for CeRu4Sbi2 consisting of the sin- 
gle valence band and the seven f states with the spin-orbit 
interaction and under the cubic crystalline field. 11 The 
former is represented by the top-most one of the valence 
bands made up of the p-orbitals on X12 clusters, and has 
the A„ symmetry which hybridizes with the f-state in the 
cluster center. The Coulomb interaction was introduced 
to all the band equally to make the use of the dynamical 
mean-field approximation possible. 

Here, we simplify the model further to investigate the 
magnetic properties of CeOs4Sbi2- The present model 
consists of a single f band (representing the J=5/2, 
ground states) and a single valence band (the above- 
mentioned A u band), hybridizing with each other. Then 
Hamiltonian is given by 

k,<7 k,tr 

+ E ( Vc lJ^ + h A + u J2 n (A 

k,(T ^ ' 

+ B ext J2<r(ni, a + 0nl <a ) (1) 

k,cr ^ ' 

where Cj CCT (/^ CT ) creates a conduction (f) electron with 
the momentum k and pseudo-spin a (representing 
doublet), and n( a =fj a f^ a with the site index i. We 
have assumed the BCC tight-binding band with near- 
est and next-nearest hopping for conduction band, e£ = 
a c ik, where = cos(k x a/2) cos(k y a/2) cos(fc z a/2) + 
a2[cos(fca;a)-r-cos(fcj,a)-r-cos(fc x a)]. The lattice constant a 
is set equal to unity in the following. a c and a-i determine 
the hopping amplitude. The f state is assumed to have 
the dispersion e£ = Ef+ctftk, namely the dispersive part 
is proportional to e£ but with a small coefficient ctf . This 
would be a reasonable choice since both X12 and R oc- 
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Fig. 1. The band structure for the non-interacting system with- 
out the external field where a c =l-0, af=0.03, a2=0.07, -Ef=0.3, 
and V=0.15. Ep represents the Fermi level. 



cupy the same BCC sites. V is the hybridization between 
conduction and f electrons. Its k-dependence is neglected. 
These choice of our model band immensely facilitates the 
calculations for the strong correlation, which will be re- 
ported in a separate paper. 12 Actually, of course, the 
f-dispersion is created through the X12 clusters, so that 
it must be self-consistently determined with e k and the 
hybridization Vk- U is the Coulomb repulsion between f 



electrons. It is well-known that e k and e 



have perfect 



nesting characters for ct2 = 0, V=0, and at half-filling. 
The last term in eq. (1) represents the Zeeman energy, 
where B ext denotes the external magnetic field along the 
z direction and (3 is coefficient of effective magnetic mo- 
ment, which corresponds to g.j( j/)=6/7x5/6=5/7 for 
Ce. 

Without the Coulomb interaction, the diagonalized en- 
ergy dispersions are obtained as = [e k + e t ± { (e k — 
e k) 2 + 4V r2 } 1 / 2 ]/2, which are the so-called hybridized 
bands, but since e£ has also a finite dispersion, the two 
bands can have finite overlap if ctf is not too small. 

Considering the empirical trend of the lattice constant 
and the energy gap of Ce-skutterudites, CeOs4Sbi2 must 
be a semimetal. 2 Our assumed band structure is shown 
in Fig. 1, which has a semi- metallic character: the top of 
lower band (T point) and the bottom of upper band (H 
point) slightly touch the Fermi level, respectively. The 
overlapping of the bands is about 0.001a c . 

Although our band structure is quite simplified, Fig. 
1 captures the low-energy structure of the band calcula- 
tion for CeOs4Sbi2 10 rather well. In the published band 
calculation, 10 one of the band made of T 8 states has a 
minimum at T point and almost touches with the valence 
band there, but it is neglected here. Also, only the single 
valence band is taken into account and the lower bands 
are completely neglected. These may affect the quanti- 
tative results but at least the qualitative feature may 
remain unchanged. The existence of T 8 band may be a 
cause of the lower temperature anomaly. 7 

In order to investigate a possibility of a magnetic 
transition and an anomalous magnetic field effect, the 
random phase approximation (RPA) is employed. Since 



Fig. 2. Static irreducible susceptibility for various momenta with- 
out the external field at T = 0.002 where x°(q) = X°i(<i) = 

X ( T \(q) = X? i (q)- 



RPA overestimates the transition temperature, we use 
the small value of the Coulomb interactions here, which 
is considered as the renormalized one at the low-energy 
regime in the rare-earth compounds. Furthermore, since 
the Coulomb interaction acts only on f electrons, ordered 
states will be mainly determined by the localized f elec- 
trons, so that only the f-f component of the correlation 
function is considered, for simplicity. It is essential for 
the following theory to calculate both the parallel and 
perpendicular elements to the external field of the static 
susceptibility, which are given by 
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where Xaa< 1S the band susceptibility, which is defined by 
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T is the temperature and N is the number of sites, fi is 
the chemical potential and u n {— (2n+ l)nT) represents 
the odd Matsubara frequency. e kcr and e k(T are the energy 
dispersion with the Zeeman terms. Note that we use the 
following relation that Xif(q) i s always equal to X|j,(q) 
within the RPA to derive \ x ■ 

We numerically calculate the above-mentioned equa- 
tions. The summations are efficiently carried out by us- 
ing the fast Fourier transform with N = 128 x 128 x 128 
points in the k summation and 4096 Matsubara frequen- 
cies in the w„ summation. Note that hopping amplitude 
a c is taken to be unity for simplicity. 

To discuss the instability to ordered states, the band 
susceptibility, x°(q) = X°t(q) = X°t(q) = X°f(q) for 
-E>cxt=0, is shown in Fig. 2. The largest value appears 
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Fig. 3. The largest susceptibilities as a function of the magnetic 
field at ' 
(X z (q)) 



Fig. 4. The T-B Phase diagram of the metal-SDW transition at 
U = 0.2. The h 
Sugawara et al' 



,.,.,„., , „,,,.,., . U = 0.2. The inset shows the experimental result obtained by 

field at T=0. 002 and (7=0.0. Solid (dashed) line represents x (q) ,2 



at H point q=(l,0, 0) in unit of 2n/a. Since the Fermi 
surfaces of electrons and holes are quite small spheres, 
the peak of x°(q) appears at q = (1,0,0) at low tem- 
peratures, which corresponds to the difference between 
the top of lower band (T point) and the bottom of upper 
band (H point) in Fig. 1. Recently, the result of the neu- 
tron scattering for CeOs4Sbi 2 at B cxt = found that the 
nesting vector corresponds to q = (1,0, 0). 13 Therefore 
our obtained result of the susceptibilities is consistent 
with the experiment at B cxt = 0. 

The basic lattice structure of CeOs4Sbi2 is BCC. For 
the single band Hubbard model at half-filling with only 
the nearest neighbor hopping, there exists the perfect 
nesting at the wave vector q = (1, 0, 0). In that case, it is 
known that the infinitesimally small U causes the metal- 
SDW transition at T=0. Although the perfect nesting is 
lost in the present model due to the hybridization and 
the next nearest neighbor hoppings, the Fermi surfaces 
are in the vicinity of the nesting. Therefore, when the 
external field is absent, the metal-SDW transition occurs 
at U c ~ 0.2 at T=0.002 within the RPA in the present 
model. 

Next, let us discuss effects due to the external mag- 
netic field. The field dependence of the largest values of 
susceptibilities is shown in Fig. 3. With increasing the 
magnetic field along the z direction, \ x is nrs t enhanced 
at small field region. As the external field is increased 
further, x x lS suppressed. The vector q which gives the 
largest values of susceptibilities does not change within 
^cxt < 0.2. In contrast to the behavior of \ x \ X z 1S 
first suppressed. This behavior is reasonable since when 
the weak external field is applied, the spin-spin correla- 
tion functions along the external field are generally sup- 
pressed. Furthermore, with increasing the external field 
further, maximum value of x z is enhanced. In this range 
of field, the vector q which gives the maximum of the sus- 
ceptibilities is gradually changed. Both of susceptibilities 
are suppressed under larger external field (not shown). 
These behaviors and the mechanisms are similar to our 
previous theory for CeRhIn 5 where more precise FLEX 



(fluctuation-exchange) method was used. 14 

In general, the response functions arc strongly affected 
by the form of Fermi surface. The present band structure 
at B cyA = is close to the perfect nesting, so that the 
change of magnetic susceptibilities becomes quite large 
for the external field. Therefore, it is possible that when 
the weak external field is applied, the nesting condition 
is rather improved because one of the spin band is shifted 
by the field and approaches the nesting condition. 

These results indicate that the presence of the ex- 
ternal field enhances instability for ordered states. In 
the 3-dimensional system, ordered states can occur at 
a finite temperature. Namely, the anomalous metal- 
insulator transition may occur in the present model with 
appropriate value of the Coulomb repulsion. Figure 4 
shows the temperature-external field phase diagram of 
metal-SDW transition. Note that the instability occurs 
not in x z (q=(l,0,0)) but in x x (q=(l, 0, 0)) (and x y ), 
which indicates the appearance of the magnetic moment 
component perpendicular to the field. It can be checked 
by an experiment in a future. In the small external field 
region, the phase transition temperature Tn is enhanced 
with increasing _B C xt , so that the slope of phase boundary 
becomes positive. Thus the anomalous phase transition 
observed in CeOs4Sbi 2 can be qualitatively reproduced. 

Note here that in order to calculate magnetic suscepti- 
bilities, we employ the RPA, which is insufficient to take 
the effect of quantum fluctuations into account. Quan- 
tum fluctuations particularly become important in the 
vicinity of quantum critical points where antifcrromag- 
netic correlations are enhanced. 14 However, correlation 
effects are essential in the weak external field region. By 
taking quantum fluctuations into account e.g. via FLEX 
approximation, 14 the metal-SDW transition temperature 
is reduced and the transition boundary as a whole shifts 
to the low temperature region. In particular, the shift 
becomes conspicuous at the weak external field region. 
However, overall structure of the phase boundary ob- 
tained here may not qualitatively changed. 

Furthermore, we comment the multi-band effects of f 
orbitals. In our model, we neglect the degeneracy of f 
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orbitals, for simplicity. However, in the published band 
calculation, one of the bands made from the 4f r§ states 
has a minimum at T point, though slightly above Ep, 
which is neglected in our model. However, we consider 
that the nesting property between the top of the valence 
band at T and the bottom of the band at H is the 
main source of the SDW transition and anomalous phase 
diagram. Therefore qualitative properties of CeOs4Sbi 2 
can be reproduced by using the present simple model. 

A similar anomalous magnetic field effect is also inves- 
tigated for the Kondo insulator, 15, 16 in which the ordered 
phase appears at finite field since the perpendicular el- 
ement of spin-spin correlations to the external field is 
enhanced. 

In summary, we investigate the anomalous phase tran- 
sition observed in CeOs4Sbi 2 . By constructing the ef- 
fective tight-binding model with the Coulomb interac- 
tion between f electrons, whose band structure is in the 
vicinity of the nesting, the metal-SDW phase boundary 
is determined within the mean field approximation. The 
nesting condition of Fermi surfaces is improved when 
the weak external field is applied, so that the slope of 
metal-SDW phase boundary becomes positive and the 
anomalous transition can be reproduced. A small mo- 
ment perpendicular to the applied magnetic field should 
be induced, which can be checked by a future experiment. 
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